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. Popular parameterizations of the freeze-out conditions in relativistic 

f^) ' heavy-ion collisions are discussed. Similarities and differences between the 

blast-wave model and the single-freeze-out model, both used recently to 
interpret the RHIC data, are outlined. A non-boost-invariant extension of 
the single-freeze-out model is proposed and applied to describe the recent 
BRAHMS data. 



PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 



1. The blast- wave model of Siemens and Rasmussen 



In 1979 Siemens and Rasmussen formulated a model describing the 
hadron production in Ne + Na F reactions at the beam energy of 800 MeV 
per nucleon pQ. The physical picture behind the model was that the fast 
hydrodynamic expansion of the produced hadronic matter leads to a sudden 
decoupling of hadrons and freezing of their momentum distributions, which 
retain their thermal character (although modified by the collective expan- 
sion effects) until the observation point. In their own words, Siemens and 
Rasmussen described the collision process as follows: "central collisions of 
heavy nuclei at kinetic energies of a few hundred MeV per nucleon produce 
fireballs of hot, dense nuclear matter; such fireballs explode, producing blast 
waves of nucleons and pions" . In this way, with Ref . , the concept of the 
blast waves of hadrons and the blast-wave model itself entered the field of 
relativistic heavy-ion collisions. 
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Although the model of Siemens and Rasmussen was motivated by an 
earlier hydrodynamic calculation by Bondorf, Garpman, and Zimanyi 
the results presented in Ref. were not obtained by solving the hydrody- 
namic equations but followed from the specific assumptions on the freeze-out 
conditions. The most important ingredient of the model was the spherically 
symmetric expansion of the shells of matter with constant radial velocity. 
With an additional assumption about the times when such shells disinte- 
grate into freely streaming hadrons (this point will be discussed in a greater 
detail in Sect. |HJ) Siemens and Rasmussen obtained the formula for the 
momentum distribution of the emitted hadrons 
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In Eq. jlj Z is a normalization factor, E = m 2 + p 2 denotes the hadron 
energy, T is the temperature of the fireball (the same for all fluid shells), 
and 7 = (1 — v 2 )- 1 / 2 is the Lorentz gamma factor with v denoting the radial 
collective velocity (radial flow). A dimensionless parameter a is defined by 
the equation 

a= —- ( 2 ) 

Small values of v (and a) correspond to small expansion rate and, as ex- 
pected, a simple Boltzmann factor is obtained from Eq. (pQ) in the limit 

v -> 0, 

S^ eXP ("^)' (3) 

The fits to the data based on the formula Q gave T = 44 MeV and 
v = 0.373. Interestingly, the value of the radial flow v turned out to be quite 
large suggesting the strong collective behavior. This was an unexpected fea- 
ture summarized by the authors with the statement: "Monte Carlo studies 
suggest that Ne + Na F system is too small for multiple collisions to be 
very important, thus, this evidence for a blast feature may be an indication 
that pion exchange is enhanced, and the effective nucleon mean free path 
shortened in dense nuclear matter" . 



2. Cooper-Frye formula 

Below we shall analyze the formal steps leading to Eq. (^P). Our starting 
point is the expression defining the momentum distribution of particles as 
the integral of the phase-space distribution function f(x,p) over the freeze- 
out hypersurface E, i.e., the renowned Cooper-Frye formula 

£ S = ^ = / d3w/(l ' I,) ' (4) 
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The three-dimensional element of the freeze-out hypersurface in Eq. (fl|l 
may be obtained from the formula 

o dx a dx^ dx^ 

dY>» = s mM ———dadf3d 1 , (5) 

where e^ a ^ is the Levi-Civita tensor and a, (3, 7 are the three independent 
coordinates introduced to parameterize the hypersurface. 

We note that for systems in local thermodynamic equilibrium we have 

E^- = j d 3 S M (x)^/ eq (^(x)^) , (6) 

where the function / e q is the equilibrium distribution function 
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Here the case e = +1 (—1) corresponds to the Fermi-Dirac (Bose-Einstein) 
statistics, and the limit e — > yields the classical (Boltzmann) statistics. 
For a static fireball one finds 

d% = (dV, 0,0,0), «„ = (1,0,0,0), (8) 

and Eq. © is reduced to the formula 

f- p = yum, o) 

where V is the volume of the system. Eq. Q agrees with Eq. © in the 
classical limit if the normalization constant Z is taken as 

V ( \i 



Z =7^^P[^)- (10) 



3. Spherically symmetric freeze-outs 

For spherically symmetric freeze-outs it is convenient to introduce the 
following parameterization of the space-time points on the freeze-out hyper- 
surface (U 

x^ = (t, x, y, z) = (i(£), r(£) sin 6 cos <f>, r(£) sin 8 sin (j), r(£) cos 9) . (11) 

The freeze-out hypersurface is completely defined if a curve, i.e., the map- 
ping £ — > (*(C)> r (C)) i n t ne t — r space is given. This curve defines the 
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Fig. 1. Left: A priori possible, different freeze-out curves in Minkowski space. The 
dotted and dashed lines describe the cases where both the space-like and time-like 
parts are present. The solid lines describe the cases where only the time-like part is 
present. Right: The (time-like) freeze-out curve assumed in the blast- wave model 
of Siemens and Rasmussen, compare Eq. I|19|) . 

(freeze-out) times when the hadrons in the shells of radius r stop to inter- 
act, see Fig. ^ The range of £ may be always restricted to the interval: 
< C < 1. The three coordinates: <p £ [0,2vr],<9 E [0,tt], and C G [0,1] play 
the role of the variables a,/3, 7 appearing in Eq. ©. Hence, the element of 
the spherically symmetric hypersurface has the form 

d 3 ^ = (r'(C), t'(C) sin 9 cos 0, t'(() sin 9 sin <f>, t'(() cos 9) r 2 (C) sin 9 d9 d<j) d(, 

(12) 

where the prime denotes the derivatives taken with respect to Besides the 
spherically symmetric hypersurface we introduce the spherically symmetric 
(hydrodynamic) flow 

utl = t(C) (1: ^(0 si 11 6* cos <ji, ?^(C) sin ^ sin 0, ^(C) cos 6 1 ) , (13) 

where 7(C) is the Lorentz factor, 7(C) = (1 — -u 2 (£)) -1 / 2 . In a similar way 
the four-momentum of a hadron is parameterized as 

= [E, p sin 9 P cos (ft p , p sin 9 P sin </> p , p cos 9 P ] , (14) 

and we find the two useful expressions: 

p. u = (E -pv(C) cos 9) 7(C), (15) 



d 3 Z-p = (Er'(() -pt'(() cos 9)r 2 (() sin 9d9d(j)d(. (16) 
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We note that the spherical symmetry allows us to restrict our considerations 
to the special case 6 P = 0. 

In the case of the Boltzmann statistics, with the help of Eqs. (@J), (fT3)) 
and (|16|) . we obtain the following form of the momentum distribution 



E 



AN 
d 3 p 



e ( E i v)/ T ' sinha dr (sinha — acosha) dt 
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(17) 

Here v , 7, r and t are functions of £, and the parameter a is defined by 
Eq. (J2J. The thermodynamic parameters T and /i may also depend on 
£• To proceed further we need to make certain assumptions about the 
^-dependence of these quantities. In particular, to obtain the model of 
Siemens and Rasmussen we assume that the thermodynamic parameters as 
well as the transverse flow velocity are constant 

T = const, n = const, v = const (7 = const, a = const). (18) 

Moreover, we should assume that the freeze-out curve in the t — r space 
satisfies the condition 



dt = v dr, t = to + v r. 
In this case we obtain the formula 

e -(Ey-l*)/T 
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jAO^dC (20) 



Equation (|2()jl coincides with Eq. 
following identification 



if we use Eq. (|TT)|) and make the 



dr 



V = -Tiri 



(21) 



Note that the quantity r max does not necessarily denote the maximum value 
of the radius of the system, see the dotted line on the left-hand-side of Fig. ^ 
An interesting and perhaps unexpected feature of the model proposed 
by Siemens and Rasmussen is the relation between the times and positions 
of the freeze-out points, see Eq. ()19|) illustrated on the right-hand-side part 
of Fig. ^ Eq. (|19|) indicates that the fluid elements which are further away 
from the center freeze-out later. Moreover, taking into account Eq. ()19|) in 
the formula (|12j) we find that the four-vector describing the hypersurface 
is parallel to the four-vector describing the flow, compare Eqs. (|12|) and 
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l|lHj) giving d 3 S^ ~ in this case. As we shall see the same features are 
assumed in the single- freeze-out model [3J E3 [7j . 

It is worth to emphasize that in the hydrodynamic approach the t — r 
freeze-out curves contain the space-like and time-like parts . The treat- 
ment of the space-like parts leads to conceptual problems since particles 
emitted from such regions of the hypersurface enter again the system and 
the hydrodynamic description of such regions (combined with the use of the 
Cooper-Frye formula) is inadeqate. Recently much work has been done to 
develop a consistent description of the freeze-out process from the space-like 
parts jSllH]. However, very often only a quantitative argument is presented 
|lf)j that the contributions from the space-like parts are small and may 
be neglected compared to the contributions from the time-like regions. The 
choice of Siemens and Rasmussen seems to have anticipated such arguments. 

4. Boost-invariant blast-wave model of Schnedermann, Sollfrank, 

and Heinz 

The model presented above is appropriate for the low-energy scattering 
processes where the two nuclei completely merge at the initial stage of the 
collision and further expansion of the system is, to large extent, isotropic. 
At higher energies such a picture is not valid anymore and, following the fa- 
mous paper by Bjorken the boost-invariant and cylindrically symmetric 
models have been introduced to describe the collisons 2 . 

The boost-invariance (symmetry with respect to the Lorentz transforma- 
tions) may be incorporated in the hydrodynamic equations, kinetic equa- 
tions, and also in the modeling of the freeze-out process. In the latter 
case, the appropriate formalism was developed by Schnedermann, Sollfrank, 
and Heinz ^2]- The ansatz for the boost-invariant, cylindrically symmetric 
freeze-out hypersurface has the form 

= (t,x,y,z) = ff(C)cosh ay, p(£) cos <j>, p(() sin cj), f(£)sinh . (22) 

Here, the parameter an is the space-time rapidity. At an = the longitudi- 
nal coordinate z is also zero and the variable f (£) coincides with the time 
coordinate t. Similarly to the spherical expansion discussed in Sect. |31 the 
boost-invariant freeze-out hypersurface is completely defined if the freeze- 
out curve £ — > (t(£), p(C)) 1S given. This curve defines the freeze-out times 

1 We use the convention that in the space-like (time-like) region the vector normal to 
the freeze-out curve is space-like (time-like). 

2 As is discussed in greater detail below, the data delivered by the BRAHMS Collabo- 
ration indicate that the systems produced at RHIC may be treated as boost-invariant 
only in the limited rapidity range — 1 < y < 1. Moreover, the assumption about the 
cylindrical symmetry is valid only for the most central data. 



7 



of the cylindrical shells with the radius p. Because of the boost-invariance 
it is enough to define this curve at z = 0, since for finite values of z the 
freeze-out points may be obtained by the Lorentz transformation. 

The volume element of the freeze-out hypersurface is obtained from 
Eq. Q, 

q .. f dp df df i dp \ 

a lu^ = I —cosh Q£|j , — cos <p, — sm0, — smh ay 1 p(QT(Qd(;da\\d<p. 

\ > > > > / 

(23) 

Similarly to Eq. (|22j) the boost-invariant four- velocity field has the structure 

u M = cosh aj_(C)cosh ay (l, tanh aj_(C) cos <p, tanh aj_(C) sin 0, tanh ay 

(24) 

We note that the longitudinal flow is simply v z = tanh ay = z/t (as in 
the one-dimensional Bjorken model), whereas the transverse flow is v r = 
tanha_i_(£)- 

With the standard parameterization of the particle four-momentum in 
terms of rapidity y and the transverse mass m±, 



p 11 = (m±cosh.y,p± cos ip, p± sin ip, mj_sinhy) , (25) 



we find 



p ■ u = m^cosh(a^)cosh(ay — y) — p^sinh(aj_) cos(0 — </?), (26) 



and 



m_LCOsh(y - ay ) ^ - p± cos(^> - (p) 



p{Qf(QdQda\\d<t>. (27) 



For the Boltzmann statistics, with [3 = 1/T, the Cooper-Frye formalism 
gives the following momentum distribution 



dN 



dyd 2 p± 



(2tt) 2 
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day / d( p(()f(() 



m_LCOsh(ay -y^jz -pj_cos(0 - ^37 



x exp 



-/?m^cosh(a^)cosh(ay — y) + /3pj_sinh(a^) cos(^ — ip) 



(28) 



The form of Eq. (|28|) shows explicitly that the distribution dN/(dyd 2 pj_) is 
independent of ?/ and in accordance with our assumptions of the boost- 
invariance and cylindrical symmetry. The integrals over ay and 4> in Eq. 
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(|28|) are analytic and lead to the Bessel functions K and I, 



dN f dp 

dyd 2 p± = 2^ m±Kl [P m ± cos H a ±)\ lo [f3p±sinh(a±)] J rf( / , (0 f (C)^ 

o 

e Pv r df 

~2tt^ P±K ° [^ m -L cosh ( a ^)] J i [/3p± sinn ( a ±)] J d ( P^(C)-^- 



o 



(29) 



In the spirit of the blast-wave model of Siemens and Rasmussen we have 
assumed here that the radial velocity is constant, v r = tanha^(C) = const, 
otherwise the Bessel functions should be kept under the integral over £. 

In order to achieve the simplest possible form of the model, the common 
practice is to neglect the second line of Eq. 1|29|). This procedure means that 
one assumes implicitly the freeze-out condition dr/d^ = (f = r = const). 
In this case the boost-invariant blast-wave model is reduced to the formula 

const m±Ki [f3m±cosh(a±)] Iq [(3p±smh(a±)] . (30) 



dyd 2 p ± 



where the constant has absorbed the factor e^ At T/)^ iax /(47r 2 ). Eq. (|3*U)) 
forms the basis of numerous phenomenological analyses of the transverse- 
momentum spectra measured at the SPS and RHIC QJ energies. 



5. Resonances 

The main drawback of the formalism outlined above is that it neglects 
the effect of the decays of hadronic resonances. Such an approach may be 
justified at lower energies but should be improved at the relativistic ener- 
gies where most of the light particles are produced in the decays of heavier 
resonance states. The expressions giving the rapidity and transverse mo- 
mentum spectra of particles originating from two- and three-body decays 
of the resonances with a specified momentum distribution were worked out 
by Sollfrank, Koch, and Heinz |141 115j . Their formulae may also be used 
to account for the feeding from the resonances in the blast-wave model, as 
already proposed in Ref. ^2]- I n other words, we wish to stress that the 
choice of the freeze-out hypersurface and of the flow profile are elements 
completely independent of the treatment of the resonances. Both are impor- 
tant with the latter being the basic ingredient in the calculation of particle 
abundances and the key to success of thermal models. 

At the SPS and RHIC energies it is important to include not only the 
decays of the most common resonances such as rj,p,uj,K* or A, but also 
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of much heavier states. Although their contributions are supressed by the 
Boltzmann factor, their number increases strongly with the mass |161 117j . 
hence their role can be easily underestimated. The effects of sequential 
decays of heavy resonances were first realized in statistical analyses of the 
ratios of hadronic abundances/multiplicities (for recent results see [T51 ITUl 
|2UJ[^) which showed that the statistical models give a very good description 
of the data, provided most of the hadrons appearing in the Particle Data 
Tables are included in the calculations. 

In order to discuss the role of the sequential decays of the resonances it 
is convenient to start with a general formalism giving the Lorentz-invariant 
phase-space density of the measured particles (22] 

ni (xi,pi) = — — = 
arp\d x\ 

I ( P2 , P1 ) [ dr 2 T 2 e-^ [ dh»W (z 2 + ^ - */ 

J E P2 J J \ m 2 

x / ( ^^-B(p N ,p N ^ 1 ) [ dT N T N e~ 

x J dT,^ (x N ) (5 (4) (x N + P ^ n ^ - XN-ij In [pn • u (x N )} . 

(31) 

Here the indices 1,2, ...,N label hadrons in one chain of the sequential de- 
cays. The first resonance is produced on the freeze-out hypersurface and 
has the label N. The final hadron has the label 1, for more details see jHlEj. 
The function B(k,q) is the probability distribution for a resonance with 
momentum k to produce a particle with momentum q in a two-body decay 

*<*•«> -iM^"**)- (32) 

The function B(k,q) satisfies the normalization condition 

(,3a B(k,q) = b, (33) 



-T]vrjv 



E q 



where b is the branching ratio for a given decay channel and p*(E*) is the 
momentum (energy) of the emitted particle in the resonance's rest frame 
(a generalization to three-body decays is straightforward and explained in 
Refs. [IH]). 

Integration of Eq. (|31[) over all space-time positions gives the formula 
for the momentum distribution 

TP dNl /V / \ 

Pl d?pl = J Xini v Cl,Pl ) = 
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B (pn,Pn-i) J dY,^ (x N ) p% f N [p N ■ u (x N )] . 

(34) 



Equation (|34|) serves as the starting point to prove that for constant 
values of the thermodynamic parameters on the freeze-out hypersurface the 
ratios in the full phase-space (4-7r) are the same as in the local fluid ele- 
ments. In this way, a connection between the measured ratios and the local 
thermodynamic parameters is obtained |23| . One may also check that for 
the boost-invariant systems it is enough to consider the ratios at any value 
of the rapidity to infer the values of the thermodynamic parameters 017]. 

The experimental RHIC data show, however, that the rapidity distri- 
butions are of Gaussian shape 3 and the thermodynamic parameters vary 
with rapidity (the measured p/p ratio depends on y), hence, the system 
created at RHIC is, strictly speaking, not boost-invariant. In this situation 
the relation between the measured ratios and thermodynamic parameters 
is not obvious. Fortunately, the RHIC data show also a rather flat rapidity 
distribution and constant ratios in the rapidity range — 1 < y < 1 |27l I28j . 
In this region (the central part of the broad Gaussian) the system to a good 
approximation may be treated as boost-invariant and the standard analysis 
of the ratios may be performed to obtain the thermodynamic parameters at 



The analysis of the ratios of hadron multiplicities measured at RHIC 
gives a typical temperature of 170 MeV. On the other hand, the analysis 
of the spectra based on Eq. (|3*U)) gives a lower temperature of about 100 - 
140 MeV. Such a situation was observed already at the SPS energies, which 
motivated the introduction of the concept of two different freeze-outs. 

Certainly, if the spectra contain important contributions from high lying 
states, the value of T obtained from the blast-wave formula fitted to the 
spectra cannot be interpreted as the temperature of the system in the precise 
thermodynamic sense. First, the contributions from the resonances (feeding 
mostly the low-momentum region) should be subtracted from the spectra of 
light hadrons, giving the insight to the properties of the primordial particles. 
Using other words, we may argue that the calculation of the ratios should 
include the same number of the resonances as the corresponding calculation 
of the spectra. 



y = 0. 



6. Single-freeze-out model 



3 This feature has revived the interest in the Landau hydrodynamic model (24II25| . see, 
for example, Ref. |26|. 
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An example of such a calculation is the single-freeze-out model formu- 
lated in Refs. EJ. In this model the decays of the resonances as well as 
the transverse flow change the spectra of the primordial particles in such 
a way that it is possible to describe well the spectra and the ratios with a 
single value of the temperature. The basic effect here is that the hadronic 
decays lead to effective cooling of the spectra. 

Similarly to the original blast-wave models discussed above, the single 
freeze-out model assumes a certain form of the freeze-out hypersurface in 
the Minkowski space. In this case it is defined by the constant value of the 
proper time 



It 2 — r 2 — r| — rf = const. (35) 
The transverse size of the system is defined by the parameter /0 maX ) 



P=\l r l + r li P</W, (36) 

and the velocity field at freeze-out is taken in the Hubble-like form 4 

M t f x y z\ 

1) T) T> T ) ■ (37) 



r t \ t t t 

The natural parameterization of the freeze-out hypersurface has the form 

t = t cosh ay cosh a±, z = rsinhaii coshaj_, 

x = t sinh a± cos <p, y = r sinh a± sin (j>, (38) 

which may be considered as the special case of the formula (|22j) . Eq. (J5J) 
leads to the following expression defining the volume element 

d£ M (x) = u^(x) t 3 sinh(a_i_)cosh(aj_) da±da\\d(f). (39) 

A very important feature of the choice ()35|) - <|37|) is that the volume 
element is proportional to the four-velocity field. This feature holds also 
in the model of Siemens and Rasmussen. In this case the treatment of the 
resonance is very much facilitated. In particular, Eq. (|34|) may be rewriten 
in the form 



1 N B (pn,Pn-i)pn -u(x n ) f N [p N -u(x N )] 



E 



dT, (x N ) pi ■ u (x N ) h\p\-u (x N )\ , (40) 



4 For a recent attempt to connect the parameterization 13511 - 1371 with hydrodynamic 
calculation see Ref. |29|. 
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where we have introduced the notation 



Pi-i ■ u (x N ) ft-i [pi_i • u (x N )\ 
d 3 p 



En 

Pi 



B(pi,pi_i)pi ■ u(x N ) fi[pi ■ u(x N )] . (41) 



In the local rest frame, the iterative procedure defined by Eq. (|41j) becomes 
a simple one-dimensional integral transform 

A-i(5) = 9^T-/ dkkfUk), (42) 
2E q p*q Jk_(q) 

where fc±(g) = m,R\E*q ±p*E q \/m\. Eqs. (JHJ) and allow us to deal 
with a very large number of decays in the very efficient way, very similar to 
that used in the calculation of the hadron abundances. 



7. Non boost-invariant single-freeze-out model 

The model described above may be generalized to the non boost-invariant 
version in the minimal way by the modification of the system boundaries. 
Introducing a dependence of the transverse size on the longitudinal coordi- 
nate z (or an), we break explicitly the assumption of the boost-invariance. 
At the same time, however, the local properties of the hypersurface and flow 
remain unchanged allowing us to treat the resonances in the same simple 
way as described in the previous Section. 

Since the measured rapidity distributions are approximately gaussian, 
it is natural to start with the gaussian ansatz for the dependence of the 
transverse size on the parameter ct\\ and restrict the region of the integration 
over aj_ to the interval 

< a± < aT ax exp[-aJ/(2A 2 )]. (43) 

The original boost-invariant version is recovered in the limit A — > oo. 
Using the values of the thermodynamic parameters obtained from the boost- 
invariant version of the model applied to the hadronic ratios measured at 
midrapidity (T = 165.6 MeV and \ib = 28.5 MeV), we are left with three 
extra parameters ( t, /?max — t sinh j^ ax , and A), which should be fitted to 
the pj_-spectra collected at different values of the rapidity. 

The result of such a fit to the available BRAHMS data on 7r + , tt~ , 
K + , and K~ production are shown in Fig. |2 The optimal values of the 
parameters found in the fit are: r = 8.33 fm, a™ 5 ™ = 0.825, and A = 3.33. 
One can see that the model reproduces the data very well in a wide range 
of the transverse-momentum and rapidity. In Fig. |31 we show the model 
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rapidity distributions compared to the data. Small discrepancies (of about 
10%) between the model and the data may be seen for the pions at y = 0. 
Note that the comparison of the rapidity distributions in Fig. |3] is done 
with a linear scale; definitely, small discrepancies may be expected for a 



BRAHMS @ 200 GeVA 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 



p T [GeV] 

Fig. 2. The single- freeze-out model fit to the transverse- momentum spectra mea- 
sured by the BRAHMS Collaboration for different values of the rapidity (2^. . The 
successive curves correspond, from top to bottom, to the rapidity bins centered 
at: y = -0.05, 0.05, 0.5, 0.7, 0.9, 1.1, 1.3, 2.2, 2.5, 3.05, 3.15, 3.25, 3.35, 3.53 EHj. 
The data in the third, fourth etc. bin are subsequently divided by factors of 2, 
which is indicated by the 2~™ label. 
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1 2 3 



y 

Fig. 3. Comparison of the measured rapidity distributions |2H] with the results of 
the single-freeze-out model. The theoretical curves were obtained by the integration 
of the spectra shown in Fig. [21 

such simplified description of the freeze-out. 

It should be emphasized that the non-boost invariant version of the 
model presented above is not capable of describing correctly the p/p ra- 
tio. In the present framework this requires an introduction of the rapidity 
dependence of the baryon chemical potential. 

8. Conclusions 

We have discussed several parameterizations of the freeze-out conditions 
in the relativistic heavy-ion collisions. We have argued that the single freeze- 
out model used to describe the RHIC data is a natural development of 
the blast-wave models worked out, among others, by Siemens, Rasmussen, 
Heinz, Schnedermann, and Sollfrank. The main advantage of the single- 
freeze-out model is that it includes all well established resonance decays, 
allowing us to treat the chemical and thermal freeze-out as essentially one 
phenomenon. In this respect, the single-freeze-out model is very similar to 
the original blast-wave model. Further similarities concern the shape of the 
freeze-out hypersurface (only the time-like parts are considered) and the 
strict use of the Cooper-Frye formula. Due to the limited space, we have 
not discussed here the variety of models where, instead of the Cooper-Fry 
formula, the so-called emission functions are introduced and modeled. An 
example of such an a approach is the Buda-Lund model j3U] . 

One of us (WF) acknowledges clarifying discussions with Jan Rafelski 
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